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1. Introduction 



Abstract. The purpose of this article is to give a complete and general answer to the recurrent 
problem in continuum mechanics of the determination of the number and the type of symmetry 
classes of an even-order tensor space. This kind of investigation was initiated for the space of 
elasticity tensors. Since then, different authors solved this problem for other kinds of physics 
such as photoelectricity, piezoelectricity, Hexoelectricity, and strain-gradient elasticity. All the 
aforementioned problems were treated by the same computational method. Although being 
effective, this method suffers the drawback not to provide general results. And, furthermore, its 

- - - complexity increases with the tensorial order. In the present contribution, we provide general 

, theorems that directly give the sought results for any even-order constitutive tensor. As an 

' — I ' illustration of this method, and for the first time, the symmetry classes of all even-order tensors 

' of Mindlin second strain-gradient elasticity are provided. 

1.1. Physical motivation. Within these last years, an increasing interest in generalized con- 
' tinuum theories [13, 11, 25, 12] has been observed. These works based on the pioneering works 

of [27, 28, 35] propose extended kinematic formulation in order to take into account size effects 
' within the continuum. The price to pay is the appearance in the constitutive relations of ten- 

■ sors of order greater than 4. These higher order objects are difficult to handle and extracting 
physical meaningful information is not straightforward. The aim of this paper is to provide gen- 
eral results concerning the type and number of anisotropic systems an even-order tensor can have. 

^ I Such results have important applications, at least, for the modeling and the nu- 

0> ■ merical implementation of non classical linear constitutive laws: 

00 ' 

• Modeling: The stakes of modeling is given a material, and a set of physical 

variables of interest to construct the more general (linear, in the present 
context) constitute law that describes the behavior of that material. An 

■ example of such a method is provided by Thionnet and Martin [32] where, 
I given a set of variables V and the material symmetry group S, they derive 

mechanical behavior laws using the Theory of Invariants and Continuum 
Thermodynamics. In such perspective our results will say, without making 
^ I any computation whether or not S is contained in the set of symmetry classes 

^ ■ of C{v, v') the space of linear applications from v £ V to v' £ V. 

- - - Numerical implementation: To implement a new linear constitutive law in a 

finite-element code, one has to known the complete set of matrices needed 
to model the associated anisotropic behavior. In that perspective, our result 
in a precious guideline as it tells you the number of the sought matrices, and 
the way to construct them. Such a way to proceed is illustrated, for example, 
in the case of 3D strain gradient elasticity in [5]. 

Constitutive tensors symmetry classes. In the field of mechanics, constitutive laws are usually 
expressed in terms of tensorial relations between the gradient of primary variables and their 
fluxes [19]. As it is well-known, this feature is not restricted to linear behaviors since tensorial 
relations appear in the tangential formulation of non-linear ones [36]. It is also known that a 
general tensorial relation can be divided into classes according to its symmetry properties. Such 
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classes are referred to as symmetry classes in the field of mechanics [15], or as isotropic classes 
(or strata) in the field of mathematical physics [1, 4]. 

In the case of second order tensors, the determination of symmetry classes is rather simple. 
According to a spectral analysis it can be concluded that any second-order symmetric tensor^ 
can either be orthotropic ([D2]), transverse isotropic ([0(2)]), or isotropic ([S0(3)]). Such kind 
of tensors are commonly used to describe, e.g., heat conduction, electric permittivity. For higher 
order tensors, the determination of the set of symmetry classes is more involved, and is mostly 
based on an approach introduced by Forte and Vianello [15] in the case of elasticity. Let us 
briefly detail this case. 

The vector space of elasticity tensors, which will be noted Ela throughout this paper, is the 
subspace of 4th-order tensors endowed with the following index symmetries: 

Minor symmetries: : Eijki = Ejiki = Ejuk 

Major symmetry: : Eijki = EkHj 
The symmetries will be denoted using the following notation: -E(jj) (ki), where (..) indicates the 
invariance under the in-parenthesis indices permutations, and ^ ^ the invariance with respect 
to the underlined blocks permutations. Index symmetries encode the physics described by the 
mathematical operator. On one hand, the minor symmetries stem from the fact that 
rigid body motions do not induce deformation (symmetry of e), and that the mate- 
rial is not subjected to volumic couple (symmetry of a). And, on the other hand, 
the major symmetry is the consequence of the existence of a free energy. An elasticity 
tensor, E, can be viewed as a symmetric linear operator on T(jj), the space of symmetric second 
order tensors. According to Forte and Vianello [15], for the classical action of S0(3), Ela is 
divided into the following 8 symmetry classes: 

[Ela] = {[1], [Z2], [D2], [D3], [D4], [0(2)], [O], [S0(3)]} 

which correspond, respectively, to the following physical classes^ : triclinic ([1]), monoclinic 
([Z2]), orthotropic ([D2]), trigonal ([D3]), tetragonal ([D4]), transverse isotropic ([0(2)]), cu- 
bic ([O]) and isotropic ([SO (3)]). The mathematical notations used above will be detailed in 
subsection 2.2. Besides this fundamental result, the interest of the Forte and Vianello paper was 
to provide a general method to determine the symmetry classes of any tensor space [4]. Since 
then, using this method, other results have been obtained: 



Property 


Tensor 


Number of classes 


Action 


Studied in 


Photoelelasticity 


T(ij)(ki) 


12 


S0(3) 


[16] 


Piezoelectricity 


T(ij)k 


15 


0(3) 


[17] 


Flexoelectricity 


T(ij)M 


12 


S0(3) 


[24] 


A set of 6-th order tensors 


Tijklmn 


14 or 17 


S0(3) 


[23] 



The limitations of the Forte- Vianello approach. The method introduced by Forte and Vianello 
is actually the more general one'^. But, in the same time, it suffers from at least two limitations: 

(1) The computation of the harmonic decomposition; 

(2) The specificity of the study for each kind of tensor. 

In its original setting, the method requires the computation of the explicit harmonic de- 
composition of the studied tensor, i.e. its decomposition into the sum of its S0(3)-irreducible 
components, also known as harmonic tensors'^. Its explicit computation, which is generally based 



Such a tensor is related to a symmetric matrix, which can be diagonahze in an orthogonal basis. The stated 
result is related to this diagonalization. 

^These symmetry classes are subgroups of the spatial group of rotations S0(3). This is due to the fact that the 
elasticity tensor is even-order. To treat odd-order tensors, the full orthogonal group 0(3) has to be considered. 

"^Some other methods can be found in the literature, such as counting the symmetry planes[10], or studying 
the SU(2)-action on Ela [8], ... but these methods are difficult to generalize to any kind of vector space. 

"^Harmonic tensors are completely symmetric and traceless, they inherit this name because of a well-known 
isomorphism in between these tensors and harmonic polynomials [6]. 
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on an algorithm introduced by Spencer [30], turns out to be intractable in practice as the ten- 
sorial order increases. But, this is not a real problem, since the only information needed is the 
number of different harmonic tensors of each order appearing in the decomposition, i.e. the 
isotypic decomposition. Based on arguments presented in [22], there exists a direct procedure 
to obtain this isotypic decomposition from the tensor index symmetries [2]. Such an approach 
has been used in [23] to obtain the symmetry classes of 6th-order tensors. 

As each kind of tensor space requires a specific study, this specificity constitutes the other 
limitation of the method. This remark has to be considered together with the observation that, 
for even-order tensors it seems that there exists, indeed, only two possibilities. Precisely, a tensor 
space has as many classes as: 

• the full symmetric tensor space: e.g. Ela is divided into 8 classes as the full symmetric 
one [15]; 

• the generic tensor space"': any other 4th-order tensor space (photoelasticity [16], flexo- 
electricity [24],...) is divided into 12 classes such as the generic one. 

The same observation can also be made for 2nd- and 6th-order tensors [23]. Understanding 
what is the general rule behind this observation would be an important result in mechanics. 
Its practical implication is the direct determination of the number and the type of symmetry 
classes for any constitutive laws no matter their orders. This result is of valuable importance 
to understand the feature of generalized continuum theories, in which higher order tensors are 
involved in constitutive laws. 

1.2. Organization of the paper. In section 2, the mains results of this paper: Theorems I, 
II and III are stated. As an application, the symmetry classes of the even-order constitutive 
tensor spaces of Mindlin second strain-gradient elasticity are determined. Results concerning 
the 6th-order coupling tensor and the 8th-order second strain-gradient tensor are given for the 
first time. Obtaining the same results with the Forte- Vianello approach would have been much 
more difficult. Other sections are dedicated to the construction of our proofs. In section 3, 
the mathematical framework used to obtain our result is introduced. Thereafter, we study 
the symmetry classes of a couple of harmonic tensors, which is the main purpose of the tool 
named clips operator. We then give the associated results for couples of SO(3)-closed subgroups 
(theorem 4.6 and table 2). Thanks to these results, and with the help of some previous works 
on that topic done by Golubistky and al. [21], we obtain in section 5 some general results 
concerning symmetry classes for general even-order tensors. In section 6 our main results are 
finally proved. The appendix is devoted to proofs and calculus of clips operations. 

2. Main results 

In this section, our main results are stated. In the first subsection, the construction of Consti- 
tutive Tensor Spaces (CTS in the following) is discussed. This construction allows us to formulate 
our main results in the next subsection. Finally, application of these results to Mindlin Second 
Strain-Gradient Elasticity (SSGE in the following) is considered. It worth noting that precise 
mathematical definitions of the symmetry classes are given in section 3. 

2.1. Construction of CTS. Linear constitutive laws are linear applications between the gra- 
dient of primary physical quantities and theirs fluxes. Each of these physical quantities (see 
table 1) are in fact related to subspaces^ of tensors spaces: theses subspaces will be called State 
Tensor Spaces (STS in the following). These STS will be the primitive notion from which the 
CTS will be constructed. In the following, C{F, G) will indicate the vector space of 
linear application from F to G. 

Let consider now two STS : Ei = and E2 = Tj, which are pth- and gth-order STS with 
possibly index symmetries. As a consequence, they belong to subspaces of (8)^R^ and i^'^M'^. A 
constitutive tensor C is a linear application between Ei and E2, i.e. an element of the space 



'The nth-order generic tensor is a nth-order tensor with no index symmetry, 
'because of some symmetries 
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Physical notions 


Mathematical object 


Mathematical space 


Gradient 


Tensor state Ti € (g)PM^ 


Tfj: tensor space with index symmetries 


Fluxes of gradient 


Tensor state T2 G (8)''M^ 


T f. tensor space with index symmetries 


Linear constitutive law 


C&C{TG,Tf) 


Tcc/:(Tg,T/) 



Table 1. Physical and mathematical links 



£(Ei, E2). This space is isomorphic, modulo the use of an euclidean metric, to Ei (8)E2. Physical 
properties leads to some index symmetries on T € Ei (8) E2; thus the vector space of such T is 
some vector subspace Tc of Ei (8) E2. 

Now, each of the spaces Ei, E2 and Ei (8) E2 has a natural 0(3) action. In this paper, we are 
concerned with cases in which p + q = 2n. In such a situation, it is known that the 0(3)-action 
on El (g) E2 reduces to the one of S0(3) [15]. We therefore have 

/:(Ei,E2) ~ El (g) E2 C TP ® T*? = TP+«=2" 

Examples of such constructions are provided in the following table: 



Property 


El 


E2 


Tensor product for CTS 


Number of classes 


Elasticity 






Symmetric 


8 


Photoelelasticity 






Standard 


12 


Flexoelectricity 


^{ij)k 




Standard 


12 


First-gradient elasticity 


^(ij)k 


T(ij)fe 


Symmetric 


17 



In the former table, two kinds of constitutive tensor spaces appeared whether they describe 

• Coupled physics : such tensors encode the coupling between two different physics, such 
as photoelasticity and flexoelectricity; 

• Proper physics : such tensors describe solely one physical phenomenon, such as classical 
and first-gradient elasticities. 

On a mathematical side this implies : 

• Coupled physics : the spaces Ei and E2 may differ, and when Ei = E2 linear applications 
are not self-adjoint; 

• Proper physics : we have Ei = E2 and linear applications are self-adjoint^. 

Therefore, the elasticity tensor is a self-adjoint linear application between the vector space of 
deformation tensors and the vector space of stress tensors. These two spaces are modeled on 
T(jj). The vector space of elasticity tensors is therefore completely determined by ^(ij) and the 
symmetric nature of the tensor product, i.e. Ela = T(jj)(8)'^T(fc;) where denotes the symmetric 
tensor product. On the side of coupling tensors, the flexoelectricity is a linear application between 
El = Tr(jj)fc, the space of deformation gradient, and E2 = the electric polarization, therefore 
Elex = T(jj)fc (g) T;. 

2.2. Symmetry classes of even order tensor spaces. Let consider an even order constitutive 
tensor space T^", it is known [22] that this space can be decomposed orthogonally^ into a full 
symmetric space and a complementary one which is isomorphic to a tensor space of order 2n — 1, 
i.e. : 

Let us introduce : 

S^": the vector space of 2nth-order completely symmetric tensors; 
G^": the vector space of 2nth-order tensors with no index symmetries^. 



This is a consequence of the assumption of the existence of a free energy. 

'The related dot product is constructed by 2n products of the canonical one. 

'Formally this space is constructed as G^" = (g)^"R^. 
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The following observation is obvious : 

and therefore, if we note 3 the operator which to a tensor space gives the set of its symmetry 
classes, we obtain: 

a(s2") c 3(t2") c a(G2") 

Symmetry group of even-order tensors are conjugate to SO(3)-closed subgroups [38, 15]. The 
classification of SO(3)-closed subgroups is a classical result that can be found in many references 
[21, 31]. These subgroups are, up to conjugacy,: 

Lemma 2.1. Every closed subgroup o/SO(3) is conjugate to precisely one group of the following 
list: 

{1, Z„, D„, r, O, I, S0(2), 0(2), S0(3)} 

Among these groups, we can distinguish: 

planar groups: : {1, Z„, D„, S0(2), 0(2)} which are 0(2)-closed subgroups; 
exceptional groups: : {T, O, X, S0(3)} which are the symmetry groups of platonician 
polyhedrons. 

Let us detail first the set of planar subgroups. If we fix a base (i; j; k) of M^, and note Q(v; 9) € 
S0(3) the rotation about v € M'^ and of angle 9 £ [0; 27r[ we have 

• 1 the identity; 

• Z„ (n > 2) the cyclic group of order n, generated by the n-fold rotation Q(k; 6* = ^), is 
the symmetry group of a chiral polygon; 

• D„ (n > 2) the dihedral group of order 2n generated by Z„ and Q(i; it), is the symmetry 
group of a regular polygon; 

• S0(2) the subgroup of rotations Q(k;6') with 9 G [0;27r[; 

• 0(2) the subgroup generated by S0(2) and Q(i;7r). 

Classes of exceptional subgroups are: T the tetrahedral group of order 12 which fixes a tetra- 
hedron, O the octahedral group of order 24 which fixes an octahedron (or a cube), and I the 
subgroup of order 60 which fixes an icosahedron (or a dodecahedron). 
In section 6, the symmetry classes of S^" and G^" are obtained: 

Lemma 2.2. The symmetry classes of S^" are : 

3(S2) = {[D2],[0(2)],[SO(3)]} 

3(8') = {[l],[Z2],[D2],[D3],[D4],[0(2)],[0],[SO(3)]} 
n > 3, a(S2") = {[1], [Z2], • • • , [Z2(„_i)], [D2], • • • , [D2„], [0(2)], [T], [O], [I], [S0(3)]} 

Lemma 2.3. The symmetry classes of G^" are : 

3{G') = {[l],[Z2],[D2],[SO(2)],[0(2)],[SO(3)]} 

3{G^) = {[1], [Z2], • • • , [Z4], [D2], ■ ■ ■ , [D4], [S0(2)], [0(2)], [T], [O], [S0(3)]} 
n>3, a(G2") = {[1],[Z2],--- ,[Z2„],[D2],--- ,[D2„],[SO(2)],[0(2)],[r],[0],[X],[SO(3)]} 

With the following cardinal properties: 



n 


1 


2 


> 3 


#a(s^") 


3 


8 


2{2n + l) 


#3(G^") 


6 


12 


4n + 5 



The symmetry classes of T are clarified by the following theorem : 

Theorem I. Let T^" be a tensor space then either ^(T^") = ^(§2") or ^(T^") = ^(G^"). 

In other terms, the number and the type of classes are the same as 

• either §2« the space of 2n-order completely symmetric tensors. In this case, the number 
of classes is minimal; 
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• or G the space of 2n-order generic tensors. In this case the number of classes is 
maximal. 

In fact, as specified by the following theorems, in most situations the number of class is indeed 
maximal: 

• For coupling tensors: 

Theorem II. Let consider T^^ the space of coupling tensors between two physics described 
respectively by two tensors vector spaces Ei and E2. If these tensor spaces are of orders greater 
or equal to 1, then 3{T^p) = 3{G^P). 

• For proper tensors: 

Theorem III. Let consider T^^, the space of tensors of a proper physics described by the tensor 
vector space E. If this tensor space is of order p > 3, and is solely defined in terms of its index 
symmetries, then J{T^p) = ^(G^p). 

Remark 2.4. Exception occurs for: 

p = 1: the space of symmetric second order tensors is obtained; 

p = 2: in the case of Tj-jj), the space of elasticity tensors is obtained. 

In each of the aforementioned situations the number of classes is minimal. There is no other 
situation where this case occurs. It should therefore be concluded that the space of elasticity 
tensors is exceptional. 

2.3. Second strain-gradient elasticity (SSGE). Application of the former theorems will 
be made on the even order tensors of SSGE. In first time the constitutive equations will be 
summed-up, then results will be stated. It worth noting, that obtaining the same results with 
the Forte- Vianello approach would have been far more complicated. 

Constitutive laws. In the second strain-gradient theory of linear elasticity [28, 14], the con- 
stitutive law gives the symmetric Cauchy stress tensor-'^'^ cr^^-' and the hyperstress tensors r^^^ 
and o;*-^-' in terms of the infinitesimal strain tensor e^^^ and its gradients r]^^^ = e^^^ (8) V and 
(X) V (8> V through the three linear relations: 

V(2) = e(4) : e(2) + m(5) .-. r/(3) + n(6) ;: ^W, 

< r(3) = M^(5) . £ + .-. r?{3) + qC^) :: k^, (2.1) 
^^(4) ^ N^(6) . ^(2) + qTH) ^(3) ^ :: 

where :,.'.,:: denote, respectively, the double, third and fourth contracted product. Above^^, 
T(ij)k, ri{ij)k = ^(ij),k, ^(ij){kl) and = ^{ij),(kl) are, respectively, the matrix 

components of cr^'^\ £^'^\ t^^\ r]^^\ uj^^^ and k^^^ relative to an orthonormal basis (i;j;k) of R^. 

And j ^(ij) (Im) , ^(ij)(/m)n> ^{ij)(kl){nin), '^ {ij)k (/m)n > 0(^ij)k(^ijn){no) and B (^ij)(^ki) (mn)(op) are the 

matrix components of the related elastic stiffness tensors. 

Symmetry classes. The symmetry classes of the elasticity tensors and of first strain-gradient 
elasticity tensors has been studied in [15] and [23]. Hence, here solely considered the spaces of 
coupling tensors N^^^ and of second strain-gradient elasticity tensors B*-^-* will be considered. 

• Let define Ces to be the space of coupling tensors between classical elasticity and second 
strain-gradient elasticity: 

Ces = {N(6) GG6|iV(,,-)(fc^)(„„)} 
A direct application of theorem 2.2 leads to the following results: 

J(Ces) = {[1], [Z2], • • • , [Ze], [D2], • • • , [De], [S0(2)], [0(2)], [T], [O], [I], [S0(3)]} 

"'^'^Exceptionally, in this subsection, tensor orders will be indicated by in-parenthesis exponents. 

^-'Xhe comma classically indicates the partial derivative with respect to spatial coordinates. The T exponent 
denotes transposition. The transposition is defined by permuting the p first indices with the q last, where p is the 
tensorial order of the image of a g-order tensor. 
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Therefore Ces is divided into 17 symmetry classes. 

• Let define Sgr to be tlie space of second strain-gradient elasticity tensors: 

Sgr = {O^*^^ G Gr^\0 {ij){kl) (mn){op) } 

A direct application of theorem 2.2 leads to the following results: 

a(Sgr) = {[1], [Z2], • • • , [Zg], [D2], • • • , [Dg], [S0(2)], [0(2)], [T], [O], [X], [S0(3)]} 
Therefore Sgr is divided into 21 symmetry classes. 

3. Mathematical framework 

In this section the mathematical framework of symmetry analysis is introduced. In the first 
two subsections the notions of symmetry group and class are introduced, meanwhile the last is 
devoted to the introduction of irreducible spaces. The presentation is rather general, and will 
be specialized to tensor spaces only at the end of the section. 

3.1. Isot ropy /symmetry groups. Let p be a representation of a compact real Lie group^^ G 

p:G^ GL{£) 

on a finite dimensional M-linear space £. This action will be noted 

y{g,x) eG X £, g-x = p{g){x.) 

For any element of £, the set of operations g in G letting this element invariant is defined as 

Sx := {5 G G I 5 • X = x} 

This set, for physicists, is called the symmetry group of x, and for mathematicians the stabilizer 
or the isotropy subgroup of x. It is worth noting that, owning to G-compactness, every isotropy 
subgroup is a closed subgroup of G. Conversely, a dual notion can be defined for G-elements. 
For any subgroup K of G, the set of i^-invariant elements in £ is defined as 

£^ ■={x££\ A; • X = X VA; € 

Such a set is referred to as a fixed point set and is a linear subspace of £. In this context we will 
note d{K) = dim £^ . It has to be observed that fixed-point sets are group inclusion reversing, 
i.e. for subgroups Ki and K2 of G, we have the following property 

KiCK2^ £^^ C £^^ 

For a given isotropy group K, the former sets are linked by the following property: 

X € ^ K C Sx 

3.2. Isotropy /symmetry classes. We aim at describing objects that have the same symmetry 
properties but may differ by their orientations in space. The first point is to define the set of all 
the positions an object can have. To that aim we consider the G-orbit of an element x of £: 

Orb(x) :={g-x\g£G}c£ 

Due to G-compactness this set is a submanifold of £. Elements of Orb(x) will be said to be G- 
related. A fundamental observation is that G-related vectors have conjugate symmetry groups. 
More precisely^^ 

Orb(x) = Orb(y) ^3geG\^^ = g^yg'^ (3.1) 
Let us define the conjugacy class of a subgroup A' C G by 

[K] = {K' C G\3g eG,K' = gKg-'} (3.2) 

In the following G will solely represent a compact real Lie group, therefore this precision will mostly be 
omitted. 

-'^''With the classical coset notation if is a subgroup of G and g £ G that is not in the subgroup H, then a 
left coset oi H in G is defined 

gH = {gh:he H} 

and symmetricaly for a right coset. 
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An isotropy class (or symmetry class) [S] is defined as the conjugacy class of an isotropy subgroup 
S. This definition implies that there exists a vector x E £" such that S = Ex and S' E [S]; 
furthermore S' = g'Sg~^ for some g & G. The notion of isotropy class is the good notion to 
define the symmetry property of an object modulo its orientation: symmetric group is related to 
a specific vector, but we deal with orbits, which are related to symmetric classes because of 3.1. 
Due to G-compactness there is only a finite number of isotropy classes [9], and we note 

3{S) :={[l];[Si];... 

the set of all isotropy classes. In the case G = SO (3) this result is known as the Hermann 
theorem [20, 2]. Elements of ^{S) are conjugate to SO(3)-closed subgroups and this collection 
was introduced in subsection 2.2. 

3.3. Irreducible spaces. For every linear subspace J- of £, we note 

g-T:= {5.x | V5 E G ; Vx E J"} 

and we say that J- is G-stable ii g ■ J- C J-" for every g G G. It is clear that, for every rep- 
resentation, the subspaces {0} and £ are always G-stable. If, for a representation p on S, the 
only G-invariant spaces are the proper ones, the representation will be said irreducible. For a 
compact Lie group, the Peter- Weyl theorem [31] ensures that every representation can be split 
into a direct sum of irreducible ones. Furthermore, in the case G = S0(3), those irreducible 
representations are explicitly known. 

There is a natural action of S0(3) on the space of M?-harmonic polynomials. If p denotes 
such a harmonic polynomial, and x E M'^, then, for every g E S0(3) we note 

g ■ p(x) = p{g~^ ■ x) 

Harmonic polynomials form a graded vector space, and to each subspace of a given degree a 
S0(3)-irreducible representation is associated. T-O' will be the vector space of harmonic polyno- 
mials of degree fc, with (Ximl-L^ = 2k+l. If we take a vector space F to be a SO(3)-representation, 
it can be decomposed into S0(3)-irreducible spaces: 

Grouping the same order irreducible spaces, the SO(3)-isotypic decomposition is obtained: 

n 
i=0 

where aj is the multiplicity of the irreducible space W in the decomposition, and n the order 
of the higher order irreducible space of the decomposition. 



Application to tensor spaces. In mechanics, F is a vector subspace of (gj^M'^. In M'^ there ex- 
ists an isomorphism, i;^, between harmonic polynomial spaces and harmonic tensor spaces[6, 15]. 
Therefore all that have been previously said in term of harmonic polynomials can be trans- 
lated in terms of harmonic tensors. A detailed discussion on this isomorphism can be found in 
[6]. Therefore U'' = ip{n''), is the space of harmonic tensors, that is the space of completely 
symmetric and traceless tensors. According to this isomorphism, any tensor space T" can be 
decomposed into S0(3)-irreducible tensors: 

n 

T" = aiM' 

1=0 
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The symmetry group of T G T*^ is the intersection of the symmetry groups of all its harmonic 
components'^ 

et = fi n ^H.. 

i=0 \j=0 / 

In the same way, rJ(T"') will be obtained as a function of the symmetry classes of the irreducible 
representations involved in the harmonic decomposition of T". The symmetry classes of S0(3)- 
irreducible representations are explicitly known [21, 18], what is unknown is how to combine 
these results to determine the symmetry classes of V (or T"). 



4. Clips operations 

The aim of this section is to construct symmetry classes of a reducible representation from 
the irreducible ones. To that aim a new class-operator, named cli'ps operator, will be defined. 
The main result of this section is given in table 2, which contains all clips operations between 
SO(3)-closed subgroups. It is worth noting that this table contains more results than strictly 
needed for the proofs of our theorems. Nevertheless, we believe that these results are interesting 
in their own and may find applications for other problems. The explicit proofs of these results 
can be found in Appendix A. 

In this section the intersection of solely two symmetry classes, is considered. Extensions 
to more general reducible representations will be treated in section 5. Let us start with the 
following lemma: 

Lemma 4.1. Let £ he a representation of a compact Lie group G that split into a direct sum of 
two G-stahle suhspaces 

£ = £i® £2 where g ■ £1 d £1 and g ■ £2 d £2^9 ^ G 

If we note by 3 the set of all isotropy classes associated to £, 3i the set of all isotropy classes 
associated to £i (i = 1,2), then [S] ^3 if and only if there exist [Si] G 3i and [II2] G ^2 such as 

s = Si nS2. 

Proof. If we take [Ei] G Ji and [S2] G J21 we know there exists two vectors xi G £1 and X2 G £2 
such that Sj = {i = 1,2). Then, let x := xi + X2. 

For every g( G Si PI S2 we have 5 ■ xi + 5 • X2 = xi + X2 = x; thus Si fl S2 C Sx. Conversely 
for every (7 G Sx we have 

5,.x = x = 5cxi+5cx2 

But, G-stable and are in direct sum, we conclude that 5 • Xj = Xj (i = 1, 2). The reverse 

inclusion is proved. 

The other implication is similar: if we take [S] G J then we have S = Sx for some x G And 
X can be decomposed into xi + X2 . The same proof as above shows that S = Sxj H Sxj • D 

Lemma 4.1 shows that the isotropy classes of a direct sum are related to intersections of 
isotropy subgroups. But as intersection of classes is meaningless, the results cannot be directly 
extended. To solve this problem, one tool named clips operator will be introduced. To that aim, 
let first consider the following lemma: 

Lemma 4.2. For every two G-classes [Sj] (i = 1,2), and for every gi, g2 in G, there exists 
g = gi^g2 in G such that 

[gi^iOi^ n 52S25^^] = [Si r\gT.2g~'^] 



In the notation H*'-' the first exponent refers to the order of the harmonic tensor, meanwhile the second 
indexes the multiplicity of in the decomposition. 
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Proof. Let g = §1^02 and 

For every 7 G S we have 7 = (7i7i5]""'^ = 92^292 ^ foi' some 7^ G Sj (z = 1, 2); then 
9il9i^ = 71 G Si and c/iT^f^ = 9729~^ ^ 9^29^^ 

Thus we have giY,g^^ C Si n §^,29^^, and conversely. As giT,g^^ is conjugate to S, we have 
proved the lemma. □ 

Definition 4.3 (Clips Operator). For every G-classes [Si] and [S2], we define the clips operator 
of [Si] and [S2], noted [Si] @ [S2], to be 

[Si] @ [S2] := {[Si n 9^29-^] for all g € G} 

which is a subset of G-classes. 

If we note 1 the identity subgroup, one can observe some immediate properties: 

Proposition 4.4. For every G-class [S] we have 

[1] @ [S] = {[1]} and [G] @ [S] = {[S]} 

Then, if we have two G-representations fi and £2, and if we note 3i the set of all isotropy classes 
of £i, the clips operator can be extended to these sets and noted 

3i®32:= U Pi]®p2] 
Sie3i, 22632 

Then, by lemma 4.1, we obtain the corollary: 

Corollary 4.5. For every two G -representations £1 and £2, if denotes the isotropy classes 
of £1 and J2 the isotropy classes of £2, then 3i ® ^2 are all the isotropy classes of £1 £2- 

Theorem 4.6. For every two SO{3)-closed subgroups Si and S2, we have 1 G [Si] ® [S2]. The 
other classes of the clips operation [Si] ® [S2] are given in the table 2 
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Table 2. Clips operations on SO(3)-subgroups 



@ 


[Zn] 


[Dn] 








[S0(2)] 


[0(2)] 


[Zm] 


[U 
















\7 1 

\u 


[ZdJ 

Z4 , [Zdz] 

[Zd] , [Drf] 












m 


[Zd,] 

[7 1 


[Z2] 

r7 1 m 1 


[Z2] 
[Z3] 

[T] 










[O] 


[Zd.] 
[Zda] 


[Z2] 

[Zda] , \Ld^ 


[Z2] 

[Zs] 

[' J 


[Z2] 

[D2] , [Z3] 
[D3] , [Z4] 
[D4] , [O] 








m 


[7 , 1 

f7j 1 


[Z2] 

[Zda] , [ZdJ 


[7 1 

r7oi 

[' J 


[7 1 

f7ol fDol 
[^3j ) [J-'Sj 

[' J 


[Z2] 

[Z3] , [D3] 
[Z5] , [D5] 
\A 






[S0(2)] 


[Zn] 


[Z2] 

[Zn] 


[Z2] 

[Zs] 


[Z2] 

[Z3] , [Z4] 


[Z2] 

[Z3],[Z5] 


[S0(2)] 




[0(2)] 


[Zd,] 


[Z2] 


[D2] 


[D2] 


[D2] 


[Z2] 


[Z2] 


[Zn] 


[Dn] 


[Z3] 


[D3] , [D4] 


[D3] , [D5] 


[S0(2)] 


[0(2)] 



Notations 

Zi := Di := 1 d2 ■= gcd{n, 2) da := gcd{n, 3) ^5 := gcd{n, 5) 

:= gcd{m, 2) dz := 2 ii d = 1, dz = 1 otherwise 

[4 if 4 I n 
"4 := < 

I 1 otherwise 

5. ISOTROPY CLASSES OF HARMONIC TENSORS 

In this section, the construction of the symmetry classes of a reducible representation from 
its irreducible components will be studied. To that aim, in subsection 5.1 the main results con- 
cerning the symmetry classes of irreducible representations are summed-up. In subsection 5.2, 
and using the results of the previous section, basic properties of reducible representations are 
obtained. These results will be used in section 6 to prove the theorems stated in section 2. From 
now on, all the results will be expressed in terms of tensor spaces. 

5.1. Isotropy classes of irreducibles. The following result was obtained by Golubitsky and 
al. [21, 18]: 

Theorem 5.1. Let S0(3) acts on M.^ . The following groups are symmetry classes ofM.^: 
(a) 1 fork> 3; 
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(b) Z„ (n > 2) for { 



n < k when k is odd 

n < — when k is even 
. - 2 

(c) Bn (n > 2) for n < k; 

(d) T fork = 3,6, 7 or k> 9; 

(e) O for k^ 1,2,3,5,7,11; 

(f ) I fork = 6, 10, 12, 15, 18 or k > 20 and k / 23, 29; 

(g) S0(2) for k odd; 

(h) 0(2) for k even; 

(i) S0(3) for any k. 

For forthcoming purposes, let us introduce the foUowing notations. For each integer k, we 
note: 

r ._ fr if rea'^ ^ __{o\iO(^3^ ^^^^ ^fxifiea^ 

^'■'^^ 1 otherwise ' ^^'^^ 1 otherwise '^^^^ 1 otherwise 

I otherwise I otherwise 

where 3^ is the set of symmetry classes of . 

5.2. Isotropy classes of direct sum. We have this obvious lemma, directly deduced from 
theorem 5.1 : 

Lemma 5.2. We have Tr^^) + ^ {[D2], [D3]} C a'^) ; V^i^^^ + ^ {[D2], [D3], [D4]} C 

; rx(fc)y^0^{[D2],--- ,[D5]}c 

We note Di(A;, n), the (n — 1) self clips operations of D''^,which is the set of isotropy classes of a 
n-uple of A:-th harmonic tensors[4], i.e. nH*^. The basic operations are, for every integers fc > 1 
and n > 2 

3{k, n) := 3'' @ 3{k, n - 1) and ^(A;, 1) := 3'^ 
On the simple example of IHI-^, the following fact can be observed 

J(2, n) := J2 ® J(2, n - 1) = 3^ @ 3^ = {[1] , [Z2] , [D2] , [0(2)] , [S0(3)]} 
This result can be generalized: 

Corollary 5.3. For every integers n>2 and k >2, the isotropy classes o/nH'^ are 

3{k, n) = 3^@3'^ = {[1], [Z2], • • • , [Z,], [D2], • • • , [Dfc], [rr(,.)], [ro(,)], [Tj^u)], P(A;)], mk)], [S0(3)]} 

Proof. From theorem 5.1 we know that [D/] G 3^ for 2 < I < k; furthermore [S0(3)] G 3^ then 
from lemma 4.4 we know that, for all integers 2 < I < k we will have (by induction), for all n > 2 
[Di] G 3{k,n). Then, when we compute 3{k,n) @ 3^ we will have [D;] @ [D;] = {[1], [Z;], [D/]} 
Neither [0(2)] and [S0(2)], with cyclic or dihedral conjugacy classes, generate other cases. The 
same occurs for clips operation of cyclic groups. Now, because of the lemma 5.2 we also see that 
no exceptional conjugacy class generates other cases. We have then prove our corollary. □ 

Corollary 5.4. For every integers 2 < 2p < 2q, we have 

3{2p, 2q) := J^P ® J^i = {[l], [Z2], ■ ■ ■ , [Z^a.(g;2p)], [D2], • • • , [D2,], 

[rT(2p) U rT(2g)], [To{2p) U ro{2q)], [rj{2p) U Ti(2q)], 
[0(2)],[SO(3)]} 

Proof. First of all, it is clear that because [S0(3)] G 3^"- {i = 1,2) we will have all [D;] for 
2<l <2q. Now we will have all [Zj] @ [S0(3)] for alll <i<p. We also have [Zj] G [D^j @ [Bj] 
for all 1 < J < 2q; this show that 

{[1],[Z2],-- - ,[Z^^(,.2p)]} ca[2p,2g] 
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Now, we can observe that clips operation of dihedral groups and [0(2)] does not generate cyclic 
groups; and lemma 5.2 shows that no other cases can be generated with exceptional subgroups. 

□ 

6. ISOTROPY CLASSES OF CONSTITUTIVE TENSORS 

6.1. The symmetry classes of even-order tensor space. Let consider the constitutive ten- 
sor space T^", it is known that this space can be decomposed orthogonally into a full symmetric 
space and a complementary one which is isomorphic to a tensor space of order 2n — 1 [22], i.e. : 

Let consider the SO(3)-isotypic decomposition of T^" 

2n 



T2" = 0«feH^ witha2n = l 

k=0 

The part related to S^" solely contains even order harmonic tensors with multiplicity one [22], 



I.e. 

n 2n—l ( „ , J, 

fcto fcto [afc -1 for fc even 

Using the clips operator, the symmetry classes of T^" can be expressed: 

j(t2") :=a(s2«)@:j(c2"-i) 

Let us first determined the symmetry classes of S^". Using the results of the previous section, 
we have: 

Lemma 6.1. 

n = l, J(S2) = {[D2],[0(2)],[SO(3)]} 

n = 2, J(S^) = {[l],[Z2],[D2],[D3],[D4],[0(2)],[0],[SO(3)]} 

V n > 3, = {[1], [Z2], . . . , [Z2(„_i)], [D2], • • • , [D2„], [0(2)], [T], [O], [Z], [S0(3)]} 

With the following cardinal properties: 

#3(S2) = 3 ; #3(S4) = 8 ; #^(§2") = 2(2n + 1) 

Proof. The case n = 1 is obtained as a direct application of the theorem 5.1 and the proposition 
4.4. For n >2, let consider the corollary 5.4 in the case ki = 2(n — 1) and k2 = 2n. 

a(2(n - 1), 2n) := {[1], [Z2], • • • , [Z2(„_i)], [D2], • • • , [D2„], [rT(2(n - 1)) U rT(2n)], 
[ro(2(n - 1)) U ro(2n)], [Ti{2{n - 1)) U rj{2n)], [0(2)], [S0(3)]} 

In the collection of planar isotropy classes, [Z2n-i] and [Z2n] are missing. It should be observed 
that the clips operation 3[2(n — l),2n] ® 3f2{n-2) never complete the sequence. 

For exceptional groups it can be observed that for any n > 3 the S0(3)-irreducible decom- 
position will contain H^. As {[7^, [O], [I]} are isotropy classes for H^, it would be the same for 
any space that contains H^. 

• Therefore, for n > 3 

3{S^n = m, [Z2], • • • , [Z2(„-i)], [D2], • • • , [D2„], [0(2)], [T], [O], [I], [S0(3)]} 

and, #3(T2") = 2(2n- 1). 

• For the case n = 2, we obtain the same result but without the class T and I and, in such 
a case, #J(Tf ) = 8. □ 

Definition 6.2. For a given S0(3) representation on the tensor space T^" {n>3), we define 

C(2n) = {[1],[Z2],-- - ,[Z2n],[D2],-- - ,[D2n],[SO(2)],[0(2)],[r],[0],[X],[SO(3)]} 
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We also define 



n = l, C(2) = {[l],[Z2],[D2],[SO(2)],[0(2)],[SO(3)]} 

n = 2, C(4) = {[1], [Z2], • • • , [Z4], [D2], • • • , [D4], [S0(2)], [0(2)], [T], [O], [S0(3)]} 



One can observe that these sets are in fact all the isotropy classes allowed by the Hermann 
theorem, and we clearly have the following cardinal properties: 

#C(2)=6 ; #(:(4) = 12 ; #C(2n, n > 3) = 4n + 5 

Definition 6.3. Let T^" be a tensor space which S0(3)-irreducible decomposition is T ~ 
0^!iQ akE^. T2" is said to be even-harmonic (EH) if a2p+i = for each < p < (n — 1). 

Lemma 6.4. Let G^" he the vector space of 2n-th order tensors with no index symmetries, G^" 
is not EH. 

Proof. For n > 1, the induced reducible SO(3)-representation on G^" = C^^^R^ is construct by 
tensorial products of the vectorial one. Such a construction implies odd-order tensors in the 
harmonic decomposition of G^". □ 

Now, the following proposition can be proved 

Theorem I. Let T^" be a tensor space, for any n > 3, if T^" is EH then ^(T^") = ^(S^"), 
otherwise ^(T^") = ^(G^"). 

Proof. Let consider the S0(3)-irreducible decomposition of T'^", this decomposition can be writ- 
ten T^" ~ S^" © C^""^. The following inclusion is always true 

a(S2") c 3(t2") C a(G2") C C{2n) 

If T^" is not EH, there exists at least one A; S N such as a2A:+i / 

3(S2") @ [S0(2)] C a(S2") @ 3(C2"-1) = 3(t2'^) 

since, from theorem 5.1, any odd-order harmonic tensor contains [S0(2)] in its symmetry classes. 
As, from lemma 6.1, dihedral groups are contained up to 2n in Di(S^"), the missing cyclic groups 
of J(§^"') are obtained by clips products with [S0(2)]. Therefore 

J(S2") ® S0(2) = C{2n) hence ^(T^'^) = C(2n) 

As G^" is not EH, rJ(G^"') = C{2n). Therefore we reach the sought conclusion, that if T^" is not 
EH a(T2") = J(G2"). 

Conversely, if T^" is EH, C^"^^ solely contains even order irreducible spaces and its leading 
harmonic spaces are, at most, of order 2(n — 1). If the leading harmonic spaces order are strictly 
less than 2{n — 1), the same analysis as for lemma. 6.1 lead to the same conclusion. Suppose 
now, that a2(„_i) > using all the previous results,we have 

a(T2") = 3(S2")@a(c2"-i) 

= (a2« @ a2("-i)) ® {3{2{n - 1), a'2(„_i)) ® J(2(n - 2), a'2(„_2))) 
= (^2" ® J{2{n - 1), 2)) = J(S2«) @ a2{n-i) 
As j2(n-i) ^Qgg contain [S0(2)] the missing classes cannot be generated, therefore J(T ) = 

a(s2"). □ 

Therefore it should be concluded that for any 2?7--order tensor space, symmetry classes are 
either the same as S'^" or the same as G'^". Let us now investigate, in the next subsection, under 
which conditions a constitutive tensor space falls into one of the identified two possibilities. 
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6.2. Construction of a constitutive tensor space. This last subsection will be devoted 
to the proof of our main result. The space of constitutive tensors is a subspace of Unear 
appUcations from Ei to E2. As seen in section 2 

£(Ei, E2) ~ El E2 C TP » T*? ~ T'^ri=p+q 

These vector spaces described the physical quantities involved in the problem understudy. We 
know, from the previous section, that any CTS has as many symmetry classes as either the 
complete symmetric tensor spaces, or the generic tensor space. Here we are interested in ob- 
taining the conditions both on Ei to E2 and on the tensor product (symmetric or not) under 
which T^*^ is even order, and therefore has a minimal number of symmetry classes. To that 
aim, distinction will be made between coupling and proper tensor spaces, in the sense previously 
defined in section 2. 

Coupling tensor spaces. We consider here two STS given by their SO (3) isotypic decomposition: 

p <? 
El = TP = and E2 = T"? = -^jW with /3p = 7^ = 1 

Now results are given by the two following lemmas: 

Lemma 6.5. //Ei / E2, p > g and ifV ® T« is EH, then V is EH and T* = EI°; 

Proof. It is sufficient to consider the tensor product of the two leading irreducible spaces, and 
to use the Clebsch-Gordan product for S0(3) [22, 2]. We obtain 

p+q 

W (^W = 

i=p-q 

Therefore p must be even, and ^ = 0, therefore T'^ = 70IHI'', and by hypothesis 70 = 1. Thus 
has to be EH. □ 

Lemma 6.6. //Ei = E2 (and then p = q) and ifT^ (g) is EH then (3i = 7^. Furthermore, if 
C{Ei) is not self-adjoint then T'P = M^; 

Proof. The demonstration is the same as the preceding proof. □ 

These results can be summed-up in the following theorem: 

Theorem II. Let consider T^p the space of coupling tensors between two physics described 
respectively by two tensor vector spaces Ei and E2 . If these tensor spaces are of orders greater 
or equal to 1, then ^(T^") = ^(G^") 

Proof. This is a direct application of the two former lemmas. □ 

Proper tensor spaces. In this case we have the following lemma: 

Lemma 6.7. //E = Ei = E2, TP (g)^ TP is EH and if C(E) is self-adjoint then 

• ifp = 2m + l then T^^+i = M^^+i; 

• ifp = 2m then T^™ = M^"^ /3oEI°; 

Proof. Because C{K) is self-adjoint the tensor product is replaced by the symmetric tensor 
product, and if TP = 0^=0 1^*^*' symmetric tensor product TP 0^ TP can be decomposed 
into a direct sum of 

PfW- and /SifSjU' (g) with i < j e {0, ■ ■ ■ ,p} 



with the following Clebsch-Gordan rule for the symmetric product : 

k 



i=0 
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Therefore, we cannot have the tensor product H* (8) HI-' for 1 < i < p — 1 and i ^ j; thus we 
deduce that 

JP = /3oIH° e BP and 0"^ = /3oIHI° (BP ®^ 11^) 2/3oHP 
Then : either p is odd, and in such case /3o = 0; either p is even and then = PqM.^ © H^. □ 

We therefore obtain the fohowing result 

Theorem III. Let consider T^^, the space of tensors of a proper physics described by the tensor 
vector space E. If this tensor space is of order p > 3, and is solely defined in terms of its index 
symmetries, then ^(T^") = ^(G^"). 

Proof. Any tensor subspace defined in terms of its index symmetries contains, as a subspace, the 
space of full symmetric tensors. Since p = 3, the harmonic decomposition of does not satisfied 
the condition of lemma 6.7. Direct application of this lemma leads to the conclusion. □ 

Remark 6.8. It can be observed, that CTS having a minimum number of classes can nevertheless 
be constructed. They consists in spaces of autoadjoint linear applications between harmonic 
spaces, which are defined both from complete symmetry under index permutations and traceless 
property. 

7. Conclusion 

In this paper the symmetry classes determination of even order tensors has been studied. 
Based on a new geometric approach, a complete and general answer to this recurrent problem 
in continuum mechanics has been given. Application of our results are direct and solve directly 
problems that would have been difficult to manage with the Forte- Vianello method. As an 
example, and for the first time, the symmetry classes of the even-order tensors involved in 
Mindlin second strain-gradient elasticity were given. To reach this goal a geometric tools, named 
the clips operator, has been introduced. Its main properties and all the products for S0(3)- 
closed-subgroups were also provided. We believe that these results may find applications in 
others context. Using the geometrical framework introduced in this paper, some extensions of 
the current method can be considered: 

• To extend this approach to odd-order tensors; 

• To take into account the coexistence of diff'erent symmetry properties for the physical 
properties of architectured multimaterials; 

These extensions will be the objects of forthcoming papers. 

Appendix A. Clips operation on SO(3)-subgroups 

In this section results concerning the clips operations on SO(3)-subgroups will be established. 
The geometric idea to study the intersection of symmetry classes relies on the symmetry de- 
termination of composite figures which symmetry groups is the intersection of two elementary 

figures. As an example let consider the rotation r = Q I k; — 1 , determining D4 n rD4r is 

amount to establish the set of transformations letting the composite figure 1 invariant. 

A.l. Parameterization of subgroups. We will define geometric elements for each S0(3) 
closed subgroup: 

• The cyclic group Z„ is characterized by the Oz axis; it will be noted Z° := Z„; 

• The same convention is retained for the dihedral group D„, i.e. := D„; 

• For the cube Co (c.f fig.2) we defined its vertex collection {Ai}i=i...s = (±1; ±1 it 1), Cq 
is O'^-invariant; 

• For the tetrahedron we consider the figure 2 and define To to be the tetrahedron AiA^AjA^, 
To is T^-invariant. 

• For the dodecahedron (c.f. fig. 3), we note "Dq the figure which vertices are^^ 

• Twelve vertices of type (±|, ib(?!)^|, O) circularly permuted; 



'With (j) the golden ratio. 
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Figure 1. Composite figure associated to D4 n rD4r* where r = Q ( k; 
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\ \ 
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\\ 








Figure 2. Cube Ci 







Figure 3. Dodecahedron Vq 



• Eight vertices of a cube with coordinates (±</'f , i*/"! , 

A. 2. Axes and subgroup classes. For every SO(3)-subgroups, we defined its g conjugate in 
the following way K'-^ = gK^g^ where the g exponent indicates the transformation, and the 
initial configuration. To proceed toward our analysis we need to introduce the following group 
decomposition [21, 18] 

Definition A.l. Let Xi, K2, • • • , K3 be subgroups of S. Then S is the direct union of the 
i^,'s if 



S = IJ if, ; h)Kir^K, = {e} Vi ^ 



J- 



i=l 



This decomposition is noted K = Ki when we have a direct union of subgroups. We 
jive some important details about geometric structure of SO(3)-subgroup. Indeed: 

• is characterized by the Oz axis, generated by k. For every rotation g S S0(3), we 
note a the axis generated by gh and note = Z^ to indicate the rotation axis. 

• DJ^ is characterized by its primary axis Oz and several secondary axis bi. Therefore 



n-l 



D°=Z°1+JZ^ 



(A.l) 



1=0 



Each bi is perpendicular to Oz and are related by the Z^ generator. is chosen such 
as one 6/ is generated by i. For every rotation g G SO (3) we define a - generated by gk 
- to be the primary axis and b - generated by gi - the secondary one, this is noted 



D 



a,b 



1 
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The subgroup can be split into a direct union of cyclic subgroups [21]. 



r° = 1+J If l+J if 

i=i j=i 



(A.2) 



where the vertex axes of the tetrahedron are noted vti and his edge axes are noted 
etj, the details of these axis appear on figure 2. Each conjugate subgroup will be 
characterized by the set of its axes {gvti, geti), g € S0(3). 
The octahedral subgroup split into 

3 4 6 

©0 = y zf . y ip i+j (A.3) 

1=1 j=i 1=1 

where vertex, edge and face axes are noted, respectively: and fcj. Detail can 

be found on figure 4. For every rotation g € S0(3), is characterized by its set of 




Figure 4. Symmetry axis of Co 

transformed axes {gfci , gecj , gvci ) . 
• The icosahedral group split into 

6 10 15 

XO = 1+J Zf ' 1+J Zf^ 1+J ' (A.4) 

i=i j=i 1=1 

where vertex, edge and face axes are noted, respectively: vdi, edj and fdj, the details 
can be found on figure 3. The vertex axes vdj are characterized by the vertices Dj for 
J = 1 • • • 10. 

A.3. Planar subgroups. 

Cyclic subgroups. First of all, we begin by this following lemma 

Lemma A.2. For every two integers m and n greater than 2, and for every two axes a and b: 

• Ifa^b then Z^ n Z^ = 1 

• If a = b then by noting d := gcd{m, n) we will have Z^ n Zj^ = Z^ 

Proof. Let g ^7J^r\ Z^, with a y^b. a and b are generated by two non-collinear eigenvectors for 
g, with eigenvalue 1. As det^ = 1 the third eigenvalue is also one, therefore g = e. Thus we 
have the first point of the lemma. If now we take, for example, a common rotation of Z^ and 
Z^, then this rotation correspond to an angle 6 = — = with r, I integers. Thus Im = rn 
and, noting m = dmi and n = dui we will have Imi = rni. As mi and ni are relatively prime, 
we deduce that 

, 1 , /, 2/7r 2aTT r. 

I = ani and then u = = — — G Zj 

n d 
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The converse inclusion is obvious, so we can conclude the lemma. □ 

A direct application of lemma A. 2 to the intersection n leads to the result: 
Lemma A. 3. For every integers n and m, we note d = gcd{n,m); we have [Z„] @ [Z^] = 

Dihedral subgroups. Let consider first the intersection T = n ZJ^. As = Z^ l+J^^^ Zg the 
following cases has to be considered: 

• When Oz = a, the intersection T = Z^n Z^ and one can apply lemma A. 2; 

• When, for some /, a = 6;, Z^ n Z'^ = 1 and one has to considered Z2 H ZJ^, which equals 
the identity as soon as m is odd. 



Lemma A. 4. For every two integers n and m, we note d := gcd{n, m) and d2{m) := gcd{m, 2); 
then we have 

[D„]®[Z^] = {[l],[Z,,(„)],[Zd]} 
Now consider the second kind of intersection 

n m 

r = n = (DO = z" 1+J z^' ) n (D^ = z^ 1+J zt ) 

1=1 1=1 

The following cases have to be considered 

• When Oz = a, and Ox = ghfi for some / : if d = 1 then T = Zg" and otherwise; 

• When Oz = a, and Ox ^ gh : T = Z°; 

• When Oz = gbi for some / : if n is even then F = Z2 and 1 otherwise. Results are the 
same when the primary axis of coincides to a secondary axis of D^. 

Lemma A. 5. For every integers n and m, we note d := gcd{n, m) and 

'2 i/d = 1 



1 otherwise 



d2{m) := gcd{m, 2) ; d2{n) := gcd{n, 2) ; dz :- 
Then we have [D„] ® [D^] = {[1], [Zrf2(„)], i'Z'diim)], ['Z'dz], [Zd], [D^] 

A. 4. Clips operations on exceptional and maximum subgroups. Here we are concern 
with the subgroups T^, O^, I^, 80(2)*^, and 0(2)'^. For these studies, we will use results con- 
cerning their proper subgroups [21]. This information is summed- up in the following diagram^^ 
[4]: 

Tetrahedral subgroup. First of all we take back the decomposition A. 2: 

7^ = zf 1 a zf^ w zf 3 w zf^ a zf' a zf^ a zf' 

We begin by the study of 7^ n Z^. As a consequence of the lemma A. 2, the primary axis of Z'^ 
must be an edge axis or a face axis of the tetrahedron. We therefore directly obain the lemma: 

Lemma A. 6. For every integer n, we note d2{n) := gcd(n,2) and d^{n) := gcd{3,n); then we 
have [Z„] © [7^ = {[1], [Z^..^^)], [Zd3(n)]} 

Now let us consider F = n . We will use primary axes and secondary axes of dihedral 
subgroup: 

We recall that the vertex axes of the tetrahedron are noted vti and his edge axes are noted etj. 
It is clear that F is a subgroup of T^. Furthermore: 

• As soon as 3 \ n a a = vti then F = Z3 is maximal; 

• When 2 I n we can find g such as a = etj, then if 6 = etj, F = D2 or F = Z2 otherwise. 

• In any case, when we only have b = etj, then F = Z2. 
Finally we can deduce the lemma: 



-'^^The arrows of figure 5 are to be understood in terms of partial order, that is inclusion of conjugates. 
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[S0(3)] 




[1] 

Figure 5. Exceptional subgroups in the poset of closed subgroup of S0(3) 

Lemma A. 7. For every integer n we note d2{n) := gcd{2,n) and d3{n) := gcd{3,n); then we 
have [D„] @[T] = m, [Z2], [Zrf3(„)], [Drf,(„)]} 

Now, for the study of F = 7^ n the arguments will be based on T subgroups as well as on 
the axes. 

• First, we can find a g such as all the axes are modified; in this case F = 1; 

• A rotation around a face or an edge axis can be found such as only this axis is leave 
fixed. Then F = Z2 or F = Z3 depending on the fixed axis. 

• If we have F D D2 then we can deduce that g carry two edge axes onto two edge axes. 
After a given permutation of axes (which leave fix T^) we can suppose that g leaves fix 
axes vti and vt2', we then conclude that g fixes also the axis vt^ and then 7"^ = T^. Thus 
we have here F = 

We deduce here the lemma: 
Lemma A.8. We have [T] @ [T] = {[1], [Z2], [Z3], [71 
Octahedral subgroup. We begin by taking back the decomposition A. 3 

3 4 6 




i=i j=i 1=1 



And, as in the case of the tetrahedron, we directly get the lemma: 
Lemma A. 9. For every integer n, we note 

{4: if 4: \ n 
1 otherwise 

Then we have [Z„] ® [O] = {[1], [Z^2(„)], [Z^.^^,,)], [Zrf^(„)] 

To study the clips operation with dihedral groups, we proceed in the same way as for the 
tetrahedron subgroup. The purpose is to examine axes of the cube and dihedral group. The 
arguments are about the same in each case. Therefore we will only detail the cases where 4 f n, 
3 I n and n is odd. 

• Either we have a = vcj, then if 6 = fcj, F = D3 or F = Z3 otherwise. 

• Or a = fcj = 0, then if 6 = fcj, F = Z2 or F = 1 otherwise. 

• Else a = ecj , then if 5 = ecj or b = fcj , F = Z2 . 
All these arguments leads us to the lemma: 
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Lemma A. 10. For every integer n, we note 

{4 if 4 \ n 
1 otrierwise 

Then we have [D„] @ [O] = {[1], [Z2], [Zd^{n)], [Z^d^n)], Pd2{n)], Pdain)], Pd^{n)]} 

Now, we take F = n and one can observe that T is necessarily a common subgroup of 
and T^, then its class must contain (in the sense of the partial order) [1], [Z2], [Z3], [D2] or 
[T]. After that: 

• There exists a rotation g around an edge axis of (that is a common face axis of the 
cube) such that only this axis is fixed; and then F = Z2; 

• There exists a rotation g around a vertex axis of (that is a common vertex axis of 
the cube) such that only this axis is fixed; and then F = Z3; 

• As soon as F D D2, as in tetrahedral case, we necessarily have F = T^. 

We conclude here the lemma: 

Lemma A. 11. We have [T] @ [O] = {[1], [Z2], [Z3], [71} 

For the study of F = 0° n O^' we will also use arguments based on subgroups. Some results 
are nevertheless more subtle: 

• First, there exists a rotation g that fix only one edge axis, and in that case F = Z2; 

• Then there exists a rotation that leaves fix only one vertex axis, and in that case F = Z3; 

• There exists also a rotation that leaves fix only one face axis, and no other axis is fixed. 
See for example figure 6(b): in that case F = Z4; 

• We can also find a rotation that leaves fix a face axis and which bring an edge axis 
onto a face axis. Indeed when we take g = Q (i; ^) we obtain F = Z^ tt) Zg = D4 (c.f. 
figure 6(c)); 

• If we take g = Q (k; f ) oQ (i; ^) we directly obtain F = D2. We can exactly compute that 
gfcs = ecg, geci = fci and gec2 = ecs and no other axes corresponds (c.f. figure 6(a)). 

• If we take g = Q(t'Ci , vr) we will find F = D3 with vci as primary axis and ec5 a secondary 
axis. 

• If F D T then, necessarily, g leaves fix the three edge axes of the tetrahedron, and then 
g will fix the cube Co; thus T = O^. 

Finally we get the lemma: 
Lemma A. 12. We have [O] @ [O] = {[1], [Z2], [D2], [Z3], [D3], [Z4], [D4], [O]} 
Icosahedral subgroup. We take the decomposition A. 4 



6 10 15 




i=i j=i 1=1 



And, as in the previous situations, we directly get the lemma: 
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Lemma A. 13. For every integer n, we note 

d2 ■= gcd{n, 2) ; ^3 := gcd{n, 3) ; ds := gcd{n, 5) 
Then we have [Z„] © [I] = {[1], [Z^^], [Za^], [Z^g]} 

Now, for the study of X° n we again use arguments about axes: 

• If a = ftj or a = vtj, then T G {Z^g, Z^^, 0^30^5} ; 

• If a = etj then F G {Za^j^di}- 

When we argue on the secondary axis of D^, we see that we can always have Z2. 
Finahy we get the lemma: 

Lemma A. 14. For every integer n, we note 

d2 '■= gcd{n, 2) ; ^3 := gcd{n, 3) ; := gcd{n, 5) 

Then we have [D„] ® [I] = {[1], [Z2], [Za,], [Z^J, [D^J, [D^g], [D^g]}} 

For the intersection n it is clear, because of the inclusion C that we can obtain 
all the classes of [T] ® [7~]. If now this intersection contains a subgroup D2, we will necessarily 
have 

D2 = Zf' W Zf' W Zf' 

where geti are the three edge axis of the tetrahedron gTo- These three axes will then have to 
correspond to three perpendicular axes of the dodecahedron. After permutation of the axes, 
which leaves fix the dodecahedron, we can suppose that these three axes are generated by the 
three vectors of the basis. But, then, the vertex axes of the tetrahedron will correspond to vertex 
axes of the embedded cube in the dodecahedron. We can then deduce that the intersection will 
be the whole T subgroup. 

Now we have to study F = Z'^ fi O^. For that, we refer to the common subgroups of [O] and 
[I]. Such subgroups can clearly be taken from the poset 5. First, it is clear that, when the cube 
related to is the embedded cube in the dodecahedron, we will have F = T. 

We also can find a rotation g such as F contains D3: indeed, g has to bring the vertex axis 
of the cube vci onto the vertex axis of the dodecahedron vd^ and the edge axis of the cube ecs 
onto the edge axis of the dodecahedron edy. With maximality argument, we can deduce that 
F = D3. We now have to examine the case of D2, Z3 and Z2: 

• When F D D2, then, after permutation of axis of the dodecahedron, we can suppose that 
g leaves fix the three axis of the basis vector. But these three axis are two order rotations 
of the dodecahedron. Thus, g will fix the cube Cq and we can deduce that = O^. We 
then have T = T- 

• We can find a rotation g around a vertex axis, for example vd^ such that F = Z3; 

• As above, we can find a rotation around an edge axis, such that F = Z2. 

We finally conclude to the formula: 

[0]@[X] = {[l],[Z2],[Z3],[D3],m} 

For the intersection F = n we will have to study the case of classes [7^, [D3], [D5], [D2], 
[Z3], [Z5] et [Z2]: 

• When F D T or F D D2, it then contains all the three two order rotations around each 
base axes, which will be three edge axes of the dodecahedron. Thus we can deduce that 
g, after permutation of these axes, leaves fix three perpendicular axes, and then g leaves 
fix X°; finally F = X; 

• There exist a rotation g around an edge axis so that F = Z2. The same argument leads 
us to Z3 and Z5. 

• If we take g to be the two order rotation around the axis ^^3 , we can compute that this 
rotation only leaves fix the axes ^^3, edg, eds and edis, and then F = D3. 

• If we take g to be the two order rotation around the face axis fdi we can also compute 
that it only leaves fix the axes fdi, ed-j, edn, edi2 and edi4, thus F = D5. 
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